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ABSTRACT 

We investigate the radial number density profile and the abundance distribution of faint satellites 
around central galaxies in the low redshift universe using the CFHT Legacy Survey. We consider 
three samples of central galaxies with magnitudes of M r = —21, —22, and —23 selected from the 
Sloan Digital Sky Survey (SDSS) group catalog of Yang et al.. The satellite distribution around 
these central galaxies is obtained by cross-correlating these galaxies with the photometric catalogue 
of the CFHT Legacy Survey. The projected radial number density of the satellites obeys a power law 
form with the best-fit logarithmic slope of —1.05, independent of both the central galaxy luminosity 
and the satellite luminosity. The projected cross correlation function between central and satellite 
galaxies exhibits a non-monotonic trend with satellite luminosity. It is most pronounced for central 
galaxies with M r = — 21, where the decreasing trend of clustering amplitude with satellite luminosity 
is reversed when satellites are fainter than central galaxies by more than 2 magnitudes. A comparison 
with the satellite luminosity functions in the Milky Way and M31 shows that the Milky Way/M31 
system has about twice as many satellites as around a typical central galaxy of similar luminosity. 
The implications for theoretical models are briefly discussed. 

Subject headings: cosmology: observations — galaxies: luminosity function — galaxies: statistics — 
Local Group 



1. INTRODUCTION 

The distribution of satellite galaxies around the central 
galaxy in a dark matter halo carries important informa- 
tion for the underlying cosmology, dark matter proper- 
ties and galaxy formation processes. It is specified by the 
abundance distribution of satellites with varying proper- 
ties (luminosity, color, star formation rate and so on), 
and the spatial distribution of the satellites. A clear pic- 
ture of how satellite galaxies are distributed in the real 
universe provides a critical test for the theoretical models 
of cosmology and galaxy formation. 

Current N-body studies have shown that the abun- 
dance of subhalos that satellites inhabit is nearly scale 
invariant, if the subhalo mass scaled by the virial 
mass of the host halo is considered (Madau et al. 2008; 
Angulo et al. 2009; Han et al. 2011). However, due to 
the complex baryonic physics involved in galaxy for- 
mation, the relation between the dark matter mass of 
a halo and its central galaxy mass is not linear (e.g. 
Yang et al. 2003; Tinker et al. 2005; Mandelbaum et al. 
2006; Wang et al. 2006; Yang et al. 2007; Zheng et al. 
2007; Guo et al. 2010; Moster et al. 2010; Wang & Jing 
2010; More et al. 2011; Li et al. 2012). Furthermore, 
when a galaxy together with its surrounding halo infalls 
into a bigger halo and becomes a satellite, physical pro- 
cesses such as tidal stripping, ram pressure and gas star- 
vation begin to take effects, which further complicates 
the galaxy-subhalo relation. Thus, the satellite distri- 
bution can not be inferred from the subhalo distribution 
nontrivially, and observations of the satellite distribution 
are expected to provide important clues to these physical 
processes. 

From rich clusters to Milky Way (MW)-sized groups, 
the luminosity of central galaxies changes by ~ 1 order 



of magnitude, while the mass of their host halos spans 
~ 3 orders of magnitude. The luminosity-halo mass rela- 
tion flattens at the massive halo end, as a result of lower 
efficiency of converting baryonic matter to light in more 
massive halos. We then expect that the number of satel- 
lites with fixed magnitude difference between the satellite 
and the central galaxy (Am = M s — M c ) is larger for 
more luminous central galaxies, unlike the nearly scale- 
independent behavior in the subhalo population. 

The Local Group satellites, due to their proximity, 
have been observed in the most detail, and the obser- 
vational data has been used to test and calibrate various 
galaxy formation models. The apparent tension between 
the large number of subhalos in ACDM simulations and 
only dozens of satellites observed in the Local Group has 
challenged theorists, making them either scrutinize var- 
ious baryonic processes that might suppress star forma- 
tion in small halos (Kravtsov et al. 2004; Koposov et al. 
2009), or investigate the possibility of alternative dark 
matter models (Zavala et al. 2009; Lovell et al. 2012). 

Despite possible underpopulation of the satellites as 
a whole in the Milky Way compared with theoretical 
subhalos, some authors have found that the existence 
of the two Magellanic Clouds (MCs) seems incompat- 
ible with the ACDM universe. Koposov et al. (2009) 
found that it is very difficult to reproduce the Magellanic 
Clouds in models that can solve the 'missing satellite' 
problem. Under the assumption of abundance matching, 
Boylan-Kolchin et al. (2010) stated that there is only a 
10% chance of finding two galaxies as bright as Mag- 
ellanic Clouds in a halo of virial mass M v - n = 10 12 M Q . 
Busha et al. (2011) reached a similar conclusion also with 
an abundance matching method. They showed that, 
even before the abundance matching is applied, the prob- 
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ability of hosting two satellites with maximum circu- 
lar velocity exceeding v max = 50kms _1 which is the 
lower bound of Magellanic Clouds (van der Marel et al. 
2002; Stanimirovic et al. 2004) is only ~ 8% for halos of 
M v ir = 1.2 x 1O 12 M in their simulation. This probabil- 
ity increases to ~ 20% for halos of M vir = 2.6 x 10 12 M Q . 
However, it is contradictory to the result of Wang et al. 
(2012), in which the number of subhalos having v max > 
60kms _1 in halos with M v i r = 2.6 x 10 12 M Q is already 
over 1. They found that the number of subhalos above 
the maximum circular velocity threshold can be under- 
estimated if the mass resolution of simulation is not high 
enough. The fact that the number of particles within 
the virial radius in Busha et al. (2011) is smaller than 
the minimum number needed to achieve convergence as 
listed in Wang et al. (2012) may account for part of the 
difference in the probability of finding two subhalos like 
Magellanic Clouds. 

On the other hand, it is doubted that the Milky 
Way, as a single case, may be representative of galax- 
ies with a comparable luminosity. Some recent works 
have been dedicated to statistical measurements of satel- 
lite abundance in large scale surveys. Guo et al. (2011) 
studied the luminosity function of satellites around iso- 
lated primaries in the Sloan Digital Sky Survey (SDSS; 
York et al. 2000), finding that the number of satel- 
lites brighter than M v = -14 in MW/M31-like galax- 
ies is lower than that in MW/M31 by a factor of 2. 
Nierenberg et al. (2012) used imaging and photometric 
redshift catalogs from the COSMOS survey to model 
the satellite distributions. Their accumulated luminos- 
ity function of satellites in MW-luminosity hosts is con- 
sistent with that in Guo et al. (2011). Robotham et al. 
(2012) searched for the Milky Way-Magellanic Clouds 
analogues in the Galaxy and Mass Assembly (GAM A) 
survey, finding that only 3.4% of MW mass galaxies have 
two companions at least as massive as the small Magel- 
lanic Cloud, and only 0.4% have an analogous MW-MCs 
system if these galaxies are all restricted to be late-type 
and star-forming. 

The spatial distribution of satellite galaxies builds 
upon that of subhalos. Many physical processes that 
affect the luminosity-halo mass relation are dependent 
on the cluster/group-centric radius. Therefore, the ra- 
dial distribution of satellite galaxies provides indispens- 
able information of various processes that have shaped 
the luminosity-halo mass relation. ACDM simulations 
have confirmed that subhalos are distributed less cen- 
trally concentrated in the host halo than the dark mat- 
ter, exhibiting a central core in the radial distribution 
of their number density (Diemand et al. 2004; Gao et al. 
2004; Ludlow et al. 2009; Gao et al. 2012). In addi- 
tion, the number density profile of subhalos has been 
shown to have a similar shape over a wide mass range. 
The radial distribution of satellite galaxies, on the other 
hand, is usually described by a power law, £(r p ) ~ 
r p ", where S(r p ) is the projected number density of 
satellites measured in observations. However, no con- 
sensus has been reached regarding the power law in- 
dex a which spans a wide range from -0.5 to -1.7 
(Madore et al. 2004; Smith et al. 2004; Sales & Lambas 
2005; Chen et al. 2006). Recent work has taken advan- 
tage of the large data sample from surveys like SDSS 



and COSMOS to study the satellite radial distribution. 
Chen et al. (2006) used the SDSS spectroscopic sample 
to study the projected radial distribution of satellites 
around isolated L* galaxies, constraining the power law 
slope to be —1.7. Watson et al. (2010) modeled the clus- 
tering of luminous red galaxies in SDSS in the framework 
of halo occupation distribution (HOD). They found that 
on projected scales 0.016 /i _1 Mpc < r p < 0.42/i _1 Mpc, 
the radial density profile of luminous red galaxies is close 
to an isothermal distribution. Nierenberg et al. (2012) 
detected satellites up to eight magnitude fainter than 
their host galaxies (with stellar mass above 10 10 ' 5 Mq) by 
subtracting the host light profile in the COSMOS survey. 
They constructed parametrized models for the satellite 
spatial distribution, number of satellites per host, and 
background galaxies, and then inferred the parameters 
using Markov Chain Monte Carlo (MCMC) method. The 
power law slope was found to be —1.1, independent of 
host stellar mass and satellite luminosities. 

Observational studies mentioned above have either 
used a spectroscopic sample, or combined a spectro- 
scopic sample for central galaxies with a deeper photo- 
metric sample to probe fainter satellites. In the latter 
case, the background must be subtracted properly due 
to the lack of redshift information. So far, several back- 
ground subtraction methods have been used. For exam- 
ple, Lares et al. (2011) applied a color cut, g — r < 1, 
to identify satellite galaxies, taking galaxies with colors 
beyond this range as background galaxies for the red- 
shift range they considered. However, this method failed 
at r p < 100 kpc. Guo et al. (2011) counted the num- 
ber of galaxies in an annuli 300 kpc < r p < 600 kpc to 
estimate the background. Nierenberg et al. (2012) used 
a similar method to build a prior in modeling the back- 
ground. Liu et al. (2011) randomized the host positions 
on the sky and searched for projected companions from 
the background. They excluded galaxies with photomet- 
ric redshift z p > 0.23 in the whole calculation to lower 
the noise, and then corrected for this photo-z loss and 
for the under-deduction of background. 

In this work, we use the halo-based group catalog of 
the SDSS constructed by Yang et al. (2007) to optimize 
the selection of central galaxies. We use a deep photo- 
metric catalog of the CFHT Legacy Survey, which is ~ 
3 magnitude deeper than the SDSS, to provide reliable 
information about the faint satellite galaxies around the 
central ones. Furthermore, we adopt a novel method of 
background subtraction that does not rely on models of 
the local environment or the photometric redshift infor- 
mation. Thus we are able to produce a reliable measure- 
ment of the satellite abundance and the radial density 
distribution for central galaxies as bright as or brighter 
than the Milky Way. 

The remaining sections are arranged as follows. Our 
methodology is described in section 2, along with the 
dataset and sample selection. Results of the radial dis- 
tributions and the abundance of satellites are presented 
in sections 3 & 4. Conclusions are given in section 5. 
Throughout this paper, we assume a ACDM cosmology 
with tt m = 0.3, fl A = 0.7, and H = lOO/ikms^Mpc" 1 
with h = 0.7. All magnitudes are given in the AB sys- 
tem. 



3 



2. DATA AND METHODS 
2.1. Data 

The data used in this work is based on the photomet- 
ric catalogue of the CFHT Legacy Survey (Gwyn 2011). 
The wide fields of the survey consist of 171 square de- 
gree pointings, which is about 145 deg 2 when the masked 
areas are excluded. The star-galaxy separation is done 
based on the Spectral Energy Distribution (SED) fit- 
ting and on the source size, closely following the method 
adopted by Coupon ct al. (2009). The survey adopts a 
5-band photometric system (u*, g' , r' , i', z'), with a limit- 
ing magnitude of r' = 25.9. To ensure a high photometric 
accuracy, we select galaxies with r' < 23.0 to be the pho- 
tometric sample from which satellites are searched for. 
In what follows, we omit the prime when referring to the 
absolute magnitude in r' band which will be written as 

We build our sample of central galaxies using the 
group catalogue constructed by Yang et al. (2007). Their 
extracted galaxy catalog is based on the NYU-VAGC 
(Blanton ct al. 2005) which is selected from the SDSS 
(York et al. 2000) data release 7 (DR7; Abazajian et al. 
2009), and supplemented with additional redshifts (e.g. 
from 2dF) to compensate for the untargeted galaxies due 
to fiber collisions. The resulting number of galaxy groups 
identified is about 4.7 x 10 5 , with spectroscopic redshifts 
0.01 < z < 0.2. After removing groups near the edges 
of the survey (see Yang et al. 2007 for details), we cross 
match the central galaxies of these groups with those in 
the overlapped area of the CFHTLS. The final matched 
sample has 3068 central galaxies, of which 82% are more 
luminous than M T — —20.5. 

2.2. Methodology 

We aim to statistically measure the number of satel- 
lite galaxies around central galaxies within a projected 
distance r p m . Given the fact that the satellites are faint 
and thus usually do not have reliable redshift informa- 
tion 1 , the number count within a projected radius around 
the central galaxies is contributed by both their satellite 
galaxies and the background galaxies 2 . Fortunately, the 
background galaxies are uncorrelated with the central 
galaxies being considered, and the contribution of back- 
ground galaxies can be statistically subtracted as shown 
in Wang et al. (2011). 

Let us consider a central galaxy at redshift z. To sub- 
tract the background contribution, we assume all galax- 
ies in the photometric survey lie at the same redshift z 
as the central galaxy. Background galaxies at redshift Zb 
with the same properties (luminosity and spectral energy 
distribution), would have the same magnitude M r , when 
these galaxies are assumed to lie at redshift z. Likewise, 
background galaxies at other redshifts may contribute to 
M r galaxies at the assumed redshift z. Since these back- 
ground galaxies with M r are not correlated statistically 
with the central one, we can then subtract the back- 
ground contribution for an ensemble of central galaxies 
by using the mean density of M r galaxies. Specifically, 

1 The photometric redshift of faint galaxies is usually very inac- 
curate for large photometry errors 

2 The foreground galaxies have the same effect as the background 
ones, so we do not distinguish them in the paper 



we generate a random sample that has the same masked 
areas and boundaries as the observed sample, and has 
random points 13 times as many as the observed sample 
has. We first count the total number of M r galaxies that 
are at a projected distance of r p away from the central 
galaxy. Next, we measure how many background galax- 
ies there would be in the random sample. Finally, we 
subtract the background contribution to obtain the sur- 
face number density of satellites S(r p ), after accounting 
for the fraction of masked area. The number of satellite 
galaxies within r p m is then 



N(< r p m ) 



Y,(r p )2irr p dr P 



(1) 



The projected correlation function w(r p ) is related to 
S(r p ) by S(r p ) = n g ii; p (r p ), where n g is the 3- 
dimentional number density of galaxies in consideration. 

We use Le Phare, a photometric redshift computing 
code (Arnouts et al. 1999; Ilbert ct al. 2006), to calcu- 
late the absolute magnitude A'/ r for each galaxy as- 
suming all satellite galaxies are at three fixed redshifts: 
z rc { = 0.06, 0.12 and 0.24. Since our central galaxy sam- 
ple spans the redshift range 0.01 < z < 0.2, we obtain 
M r at redshift z by 

M T {z) = M r (z ie {)+(DM(z Ie i)-DM(z))+(kcor(z Ie {)-kcor(z)), 

where DM is distance modulus. We compute the differ- 
ence in the /e-correction, kcor(z re {) — kcor(z), by linearly 
extrapolating the fc-correction between the two adjacent 
reference redshifts. 

Galaxies are distributed in the form of clustering, with 
more luminous galaxies being more strongly clustered. 
Galaxies of characteristic luminosity L* show a correla- 
tion length of ~ 5 ft. _1 Mpc (see, e.g. Zehavi et al. 2005; 
Li ct al. 2006), which is much larger than the virial radius 
of their host dark matter halos. Therefore, the number 
excess of satellites around central galaxies after deduct- 
ing a uniform background comprises two components: 
one from satellites within the virial radius, and one from 
the correlated nearby structures outside of the virial ra- 
dius along the line of sight. The projected correlation 
function thus consists of two parts, 



w p (r p ) = w pl (r p ) + w p2 (r p ), 



where 



and 



w p i(r p )=2 / rdr£(r)(r - r p ) 



2 _ - 2N.-1/2 



w p2 (r p ) = 2 rdr£{r)(r - r p ) 



2 _ r 2\-l/2 



(3) 



(4) 



(5) 



Here £(r) is the real space correlation function 
(Davis & Peebles 1983). If £(r) follows the power law 
form £(r) = (r/7' )~ 7 , then the overall projected correla- 
tion function becomes 



7-1 



/r 



(6) 



The contribution of real satellites (galaxies within r v ; r ) 
to the surface density at the projected distance r p can 
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thus be obtained by 

where Wpl ^ Tp ^ is a decreasing function of r p because of 

the projection effect. Therefore, the total number of real 
satellites within the projected distance r p m would be 

JV1(< r p m ) = T P £(r p )^^27rr p dr p . (8) 

JO w P\ r p) 

According to equation (8), we need to know the ratio 
of Wpl ( Tp J to obtain the number of real satellites within 

the projected distance r p m . The ratio of -^j^y, which 

can be calculated with equations (4) and (6), is inde- 
pendent of the correlation length r , and depends only 
on the power law index 7. The value of 7 can be in- 
ferred from the projected radial distribution profile of 
satellites. The details arc described in section 3.1. There 
we find a power law index of —2.05 for central galaxies in 
the three magnitude bins we consider (M r c = —21, —22, 
and —23), irrespective of satellite luminosities. With 
7 = —2.05, for central galaxies with M t = —21,-22, 

and —23, the percentage j^t^t^t of real satellites within 

the projected distance r p m = 0.3 Mpc is 66%, 79% and 
88% respectively. The number of real satellites within 
the projected virial radius iVj(< r v ; r ) is equivalent to the 
number of satellites within the virial radius, which is 66% 

Of N{< 7Vir). 

3. RADIAL DISTRIBUTIONS OF SATELLITES 

In this section we study the projected radial distribu- 
tion of satellites. We first consider the projected number 
density profile, and we fit it with a power law to infer 
the logarithmic slope 7 of £(r) in order to feed the the 
method elaborated in section 2.2. We then calculate the 
projected two-point correlation function w p (r p ) to inves- 
tigate how the correlation function changes with satellite 
luminosity. 

3.1. Projected surface density profile 

We consider three central samples with magnitudes in 
bins centered at M T ° = —21,-22, and —23 with an in- 
terval of 1 magnitude. In the SDSS DR7 group cata- 
logue, host halo masses are estimated using the rank- 
ing of the groups according to their characteristic group 
luminosity or stellar mass. Yang et al. (2007) showed 
that, the characteristic stellar mass is more tightly cor- 
related with the halo mass than the characteristic group 
luminosity. Therefore, we choose the halo mass esti- 
mated from the stellar mass to calculate the virial ra- 
dius. In the faintest luminosity bin with M T C = —21, 
the halo masses are about 90% complete. Since this lu- 
minosity is close to that of the Milky Way, we consider 
Milky Way-sized host halos: 1.0 xlO 12 M < M ha i <3.0 
xl0 12 M Q (e.g. Wilkinson & Evans 1999; Sakamoto et al. 
2003; Dehnen et al. 2006; Xue et al. 2008), and draw the 
virial radius from a cosmological N-body simulation. The 
simulation is the same as that used in Jiang et al. (2008), 
and the details can be found there. The resulting virial 



radius is 0.3 Mpc. Those of the two more luminous bins 
are 0.43 Mpc and 0.73 Mpc respectively. 

To see how the distribution changes with the lumi- 
nosities of both the centrals and the satellites, we first 
consider three subsamples of satellites for each central 
galaxy sample: M r s < M r c + AM, where AM = 2,3,4. 
We count satellites from 0.01 Mpc to 1 Mpc away from 
their central galaxies. We scale the surface number den- 
sity by the mean value within the virial radius, and the 
distance by the virial radius. These radial distributions 
are displayed in the logarithmic space in Figure 1. The 
errors are estimated by dividing each sample of central 
galaxies in the specific magnitude bin into eight subsam- 
ples. Central galaxies are randomly grouped to form sub- 
samples. 

For the most luminous sample of central galaxies, the 
surface density drops at log(r p /r v ; r ) < —1.5 for all the 
three satellite samples. It corresponds to a projected 
distance of < 16 kpc, which is even smaller than the ra- 
dius of galaxies as luminous as M r c = —23. Therefore, 
satellites could be cannibalized by the central one if their 
separations become comparable to this scale, so that the 
projected surface density is lowered. For central galaxies 
of M r c = —21, the surface density rises once satellites 
are beyond the virial radius. These results are consistent 
with the outer structure of dark matter halos found in 
pure N-body simulations (Prada ct al. 2006, e.g.), where 
they showed the overdensity of dark matter around halos 
of mass < 5 x 10 12 Mq actually rises with the distance at 
« 2 — 5r vir . The fact can be understood, because groups 
with mass < 5 x 10 12 M Q are located along filaments con- 
necting more massive halos, unlike massive clusters (or 
halos) which are formed at density peaks. It is amaz- 
ing to see that the rising amount found for the groups of 
M T C = — 21atl — 3 r v ir is in good agreement with that 
found by Prada et al. (2006), although the error bar is 
large at r p > r vir . 

Within the virial radius, the density profiles of the 
satellites are well approximated as linear relations in 
the logarithmic space. Therefore, as in previous works 
(Chen et al. 2006; Chen 2008; Nierenberg et al. 2012), 

we fit the radial distribution with a power law, s ^. p ^ ^ = 
A(-p p -) a , using data in the range —1.5 < log(?' p /r v i r ) < 
0, where A = (a + 2)/2. The best-fit parameter a is 
listed in Table 1. All the slopes are close to —1.0 with 
a small scatter. They are independent of both the lumi- 
nosity of central galaxies and that of satellite galaxies. 
The median value of a is — 1.05 ± 0.08, very close to the 
isothermal value. The blue dotted line in each panel of 
Figure 1 shows the power law function with a = —1.05. 

It is interesting to compare the satellite distribu- 
tion with the subhalo distribution or the dark mat- 
ter distribution in the dark matter halos. The den- 
sity profile in the dark matter distribution in halos has 
been studied in many works (e.g. Navarro et al. 1997; 
Ghigna ct al. 2000; Jing & Suto 2000; Klypin et al. 2001; 
Navarro ct al. 2004). Among them, the most widely used 
is the NFW profile (Navarro ct al. 1997). In an NFW 
profile, the logarithmic slope changes slowly from —3 in 
the outer region to approaching —1 in the inner region 
of the halo. In contrast, subhalos do not follow the same 
distribution, but have a much shallower density distribu- 
tion than the dark matter density profile (Dicmand ct al. 
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2004; Gao et al. 2004, 2012). Interestingly, these authors 
have shown the subhalo radial distributions are nearly 
independent of the subhalo mass when the simulation 
resolution is properly taken into account. For compar- 
ison, we plot typical dark matter and subhalo density 
profiles in Figure 1, taken from Han et al. (2011) af- 
ter the projection. Around the radius close to the virial 
radius, the logarithmic slope of subhalos is about — 1, 
very close to the observed value of satellites, but much 
flatter than the dark matter density profile which has a 
slope —3 in 3-dimentional space. This is expected be- 
cause the subhalos around the virial radius have suffered 
from relatively weak tidal stripping of the host halo, and 
the subhalo mass has a close correspondence to the lu- 
minosity of satellites. In the inner part of the halo, the 
subhalos suffer from a stronger tidal stripping, with the 
strength increasing as the distance from the halo center 
decreases. This leads to a steeper radial distribution of 
the satellite galaxies than that of subhalos. Although the 
satellite distribution can be qualitatively explained, the 
unique non-trivial features of the radial distribution of 
satellites as shown in Figure 1 should provide a quanti- 
tative test for theories of galaxy formation, especially on 
the processes of tidal stripping, dynamical friction, and 
gas ram pressure. 

Next, we consider another two subsamplcs of fainter 
satellites: M r s < M r c + 5 and M r s < M r c + 6. Figure 
2 compares radial distributions of these two subsamples 
with the distribution function obtained above. The ra- 
dial distributions for fainter satellites still conform to the 
power law form at large radii. It begins to deviate from 
the power law at small radii. It drops sharply for all 
central samples in the subsample with M r s < M r c + 6. 
The deviation appears at r p ~ 0.2r v j r for the two lumi- 
nous central samples. For the faintest central sample, 
the deviation appears even at a larger distance. 

There are two possible reasons why the density of faint 
satellites drops at small scales. The first one is observa- 
tional. Faint satellites close to a much brighter central 
galaxy may be missed from the photometric catalogue. 
The second reason is physical. In the central region, 
subhalos that host satellites are heavily stripped by the 
tidal force. This leaves the satellites subject to the tidal 
stripping. The fainter the satellite is, the more it may 
be stripped. Nierenberg et al. (2012) have examined the 
faint satellite distribution around bright centrals by care- 
fully removing the bright central background with the 
quality photometric image of the COSMOS field. They 
found that the logarithmic slope is —1.1 even at the dis- 
tance close to the bright central galaxy. With their find- 
ings, we think that the drop of the faint satellite density 
around the central ones is more like caused by the obser- 
vation or the data reduction. One needs to take this into 
account when computing the luminosity function and 
the abundance of satellites in observations. Considering 
the fact that the surface density of the faintest satellites 
starts to drop at r p /r V j r ~ 0.2 for the M r = —21 central 
sample, we estimate that both observational quantities 
(the luminosity function and the abundance) could be 
underestimated by 20 percent. This effect is small for 
other cases of satellites and central galaxies. 

3.2. Projected two-point cross correlation function 



Based on the power law radial distribution function 
obtained above, we obtain a smoothed surface number 
density £(r p ) of satellite galaxies. Then we adopt a 
method similar to Wang et al. (2011) to calculate the 
projected two-point cross correlation function w p (r p ) of 
the central galaxies and the satellites, based on the rela- 
tion £(r p ) = w p (r p )n g . We estimate n g from the r'-band 
luminosity function, which is approximated by that of the 
SDSS r o f -band given in Blanton et al. (2003). 

Figure 3 shows the projected cross correlation func- 
tions between central galaxies in three magnitude bins 
and satellites in six relative magnitude bins. The loga- 
rithmic slope is the same for all cross correlation func- 
tions, as expected from the number density profile. The 
amplitude of w p varies both with the luminosity of cen- 
tral galaxies and satellite galaxies. A conspicuous feature 
displayed in the plots is that, the clustering amplitudes 
do not exhibit a monotonic dependence on satellite lu- 
minosity. This is most pronounced for the faintest cen- 
tral magnitude bin. The amplitude first decreases as 
the satellite luminosity is lowered to M r s — M r c = 2, 
and then it increases as the satellite luminosity continues 
to be lowered. The turning point is relatively at larger 
AM = M/ — M r c in the most luminous central sample. 
This happens at M r s - M r c = 4. The relatively close 
amplitudes between different satellite bins for central 
galaxies with M r c = —22, —23 are generally consistent 
with those found in Wang et al. (2011). Their sample 
is not deep enough, however, so that the non-monotonic 
trend cannot be detected in their work. A similar be- 
havior was found in Li et al. (2006) and Zehavi et al. 
(2011). Li et al. (2006) found that, the projected auto- 
correlation function for the red galaxies at r p = 0.2/i _1 
Mpc reaches the lowest value when the galaxy luminosity 
is around —20.5+5 log 10(h). Zehavi et al. (2011) showed 
with their sample of red galaxies that, the decreasing 
trend of projected auto-correlation functions with lumi- 
nosities is reversed at r p <~ 2/i -1 Mpc when the galaxy 
luminosity is fainter than —20 + 51ogl0(/i). Since the 
cross correlation is an approximation of the geometric 
mean of the two auto-correlations (Szapudi et al. 1992), 
the rebound behavior in the auto-correlation functions 
indicates a similar feature in the cross correlations. The 
higher clustering amplitudes of faint satellite galaxies (es- 
pecially those with M r s — M x ° ~ 5.5) around their cen- 
trals suggests a higher fraction of these galaxies exists as 
satellites in halos being considered than the more lumi- 
nous galaxies. 

4. LUMINOSITY FUNCTIONS OF SATELLITES 

Figure 4 shows the luminosity functions of satel- 
lites for central galaxies in the three magnitude bins: 
M r c = -21,-22, and -23. The number of satellites 
per magnitude for each central galaxy is plotted against 
the r' band magnitude difference between satellites and 
their central galaxies, AM = Al/ — M T C . As in section 
3.1, we estimate the errors from eight randomly grouped 
subsamples. Satellites are counted within the virial ra- 
dius of their host halo, and the number count plotted in 
this figure and the following figures are iVj (see Equation 
8). 

With Figure 4, we compare the population of satellites 
with that of subhalos to find their connections. Using the 
scaled subhalo mass which is the subhalo mass divided 
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by their host halo mass, M su b/Mh os t, studies based on N- 
body simulations have found that the subhalo mass func- 
tion dN/dlog^^ oc (i|f;) /3 with $ « -0.9, depend- 
ing on the host halo mass very weakly (at most Mhost j 
Angulo et al. 2009; Han et al. 2011). If the mass-to-light 
ratio is the same for all (sub)halos, we can easily deduce 
that the number of satellites per central galaxy when 
expressed function of AM as in Figure 4 should 
also depend little on the luminosity of the centrals, but 
increase linearly with AM. In contrast, Figure 4 demon- 
strates a different picture. The major contribution of the 
difference between the three magnitude bins comes from 
a varying mass-to-light ratio with halo mass. We know 
that, the mass-to-light ratio at z ~ is the lowest at 
Mhaio ~ 10 12 M Q (Yang et al. 2003; Zheng et al. 2007; 
Zehavi et al. 2011), which marks the halo mass with the 
highest star formation efficiency. Above this mass, the 
mass-to-light ratio rises and the light increases more and 
more slowly with mass. Therefore, the number of satel- 
lites at fixed AM is different for central galaxies with 
different luminosities. It depends on the difference in 
the mass-to- light ratio between the central and the satel- 
lite galaxies. This difference in mass-to-light ratio does 
not solely come from the difference in star formation ef- 
ficiency, but also from the fact that satellite halos are 
stripped of their mass after being accreted to host halos, 
leading to a lower mass-to-light ratio compared to that 
at accretion time. 

To compare with a previous work by Guo et al. (2011), 
we also calculate the luminosity functions within 0.3 Mpc 
of central galaxies (solid lines in Figure 5). Results from 
Guo et al. (2011) are shown with dotted lines. We see 
that, our results are higher than theirs in all luminosity 
bins. Although the r' filter we use is slightly different 
from the SDSS r filter, according to the color transfor- 
mations with the SDSS filters (Rcgnault et al. 2009), we 
find a tiny difference of 0.01 mag between r' and tsdss 
for central galaxies we consider here. One reason for the 
discrepancy is the different method adopted in the sub- 
traction of the background galaxies. Guo et al. (2011) 
take galaxies that are 0.3 Mpc to 0.6 Mpc away from 
their central galaxies as a proxy for background galax- 
ies. However, this distance range is well within the virial 
radius for galaxies of M r c = —23, and partly within the 
virial radius for galaxies of M r c = —22, according to the 
virial radius we calculated above. Therefore, the back- 
ground may be over-subtracted for these two luminosity 
bins. For central galaxies of M r c = —21, although the 
radial range of 0.3-0.6 Mpc already falls out of the virial 
radius, it is however not guaranteed that the background 
galaxies are properly subtracted in this way. In fact, as 
shown by Figure 1 and the N-body results of Prada et al. 
(2006) , the nearby structures outside the virial radius of 
the host halo of a Milky Way-like galaxy have elevated 
the surrounding galaxy density. We also note that in 
their work, the correlated galaxies outside the virial ra- 
dius along the line of sight are not explicitly subtracted. 

Next, we study the satellite abundance in MW/M31- 
like galaxies. The V-band magnitude (vega system) of 
the Milky Way and M31 is -20.9 and -21.2 respectively 
(van den Bergh 2000), which gives a mean magnitude 
of -21.05 for the MW/M31 system. For galaxies with 
M v = —20.9 ± 0.5 (the conversion from vega magnitude 



to AB magnitude is ignored, which is only 0.02), the 
mean magnitude in the r' band is M r = —21.3. Liu et al. 
(2011) discussed that M v = -20.9 ± 0.5 corresponds to 
M r 01 = -21.2, in which M^- 1 is the SDSS r band mag- 
nitude which is fc-corrected to z = 0.1. Considering the 
fc-correction from z = 0.1 to z = 0.0 is of order 0.1, the 
two values are consistent. In addition, for galaxies with 
M v = -21.05 ± 0.5, the mean M r = -21.4. We consider 
both samples, one with M v = —20.9 ± 0.5, and the other 
with M v = —21.05 ± 0.5. We convert the r' band to the 
V band magnitude with the V — r' color. Then we search 
for the centrals in the M v interval from all the matched 
central galaxies. The difference of the two samples are 
very small, as shown in Figure 6. We can see from Fig- 
ure 6 that, the luminosity function of bright satellites 
in the Milky Way and M31 (red points, Koposov et al. 
2008) lies above that of our MW/M31-like galaxies by a 
factor of 2. A comparison with the results in Guo et al. 
(2011) is also given in this figure. It shows a reasonable 
agreement except that the amplitude around AM = 3 is 
a bit lower in their results, despite the different methods 
used. We also note that the satellite luminosity functions 
of the two Milky Way like galaxy samples (M v = —21.05 
and My = -20.9) are close to that of M r c = -21 galax- 
ies in Figure 5 in our analysis as expected, but they are 
not apparently in Guo et al. (2011) (comparing the green 
dotted line in Figure 5 with the blue line Figure 6). 

5. CONCLUSIONS 

In this work, we have investigated the distributions of 
satellites around their central galaxies using the photo- 
metric catalogue of the CFHT Legacy Survey. We use 
the SDSS/DR7 group catalog of Yang et al. (2007) to 
identify the central galaxies. The method we have used 
can count in all the candidate satellites while subtract- 
ing the background galaxies accurately. We have focused 
on three samples of central galaxies with magnitudes of 
M r c = —21,-22, and —23. Our main results can be 
summarized as follows. 

• Over a wide radial range — 1.5 < log(?' p /r vir ) < 0, 
the distribution of projected radial number den- 
sity obeys a power law form with the best-fit log- 
arithmic slope of —1.05, indicating that the 3- 
dimensional number density follows approximately 
an isothermal distribution. The radial number den- 
sity profile, if scaled by the mean number density 
within the virial radius, is independent of the lu- 
minosity of central galaxies and satellite galaxies. 

• The projected cross correlation functions between 
central galaxies and their satellites exhibit a non- 
monotonic trend with satellite luminosity. As the 
satellite luminosity decreases, the amplitude of the 
cross-correlation function first decreases, before it 
increases after the satellite luminosity crosses a 
turning point. This is most pronounced for central 
galaxies with M r c = —21, where the decreasing 
trend of the clustering amplitude as the satellite 
luminosity decreases is reversed when satellites are 
fainter than central galaxies by more than 2 mag- 
nitudes. 

• The Milky Way/M31 system has about twice as 
many satellites with magnitude difference AM = 



M T S — M T C < 6 as we have obtained statistically 
for galaxies of the same luminosity. This over- 
abundance shows that the MW/M31 is atypical of 
central galaxies with the same luminosity in the 
distribution of satellites. 
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Fig. 1. — The surface density of satellites as a function of the distance to the central galaxies. The surface density is scaled by the 
mean surface density within the virial radius, and the distance is scaled by the virial radius. From left to right: M r " < M r c + 2, M r a < 
M r c + 3, M r s < M r c + 4. The blue dotted line is a fitting line with a slope of -1.05. The radial distributions for dark matter (cyan solid) 
and subhalos (purple solid) are also shown in the third panel. For clarity, error bars are plotted only for four subsamples, and those for 
other subsamples are comparable. 
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Fig. 2.— Similar to Figure 1, but for satellites with M r s < M r c + 5 and M r s < M r c + 6. 
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Fig. 3. — Projected correlation functions for central galaxies with M r c = 



—21, —22, —23 from left to right, extending to virial radii. 




Fig. 4. — Luminosity functions for central galaxies in the three magnitude bins within virial radii. 




Fig. 5. — Similar to Figure 4, but satellites are counted within projected distance 0.3 Mpc. In comparison, results of Guo et al. (2011) 
are shown with dotted lines. 
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Fig. 6. — Luminosity functions for Milky Way-like central galaxies (M v = —20.9 ± 0.5) (black dotted) and Milky Way/M31-like central 
galaxies (M v = —21.05 ± 0.5) (black solid). The luminosity function of the Milky Way and M31 is plotted with red points (Koposov ct al. 
2008). The results of Guo et al. (2011) are shown with the blue solid line. 



